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^::j- ' The boost-invariant tunneling of particles along the hyperbolas of constant invari- 

^N ■ ant time r = \/t'^ — z^ is included in the transport equations describing formation of 

the quark-gluon plasma in strong color fields. The non-trivial solutions of the trans- 

^ \ port equations exist if the boost-invariant distance between the tunneling particles, 

^^ ■ measured in the quasirapidity space t] = l/2ln{{t-\- z) / {t — z)), is confined to a finite 

pg . interval Arj. For realistic values of Arj the solutions of the transport equations show 

C^ I similar characteristics to those found in the standard approach, where the tunneling 

takes place at constant time t. In the limit Arj -^ cxd, the initial color fields decay 

instantaneously. 
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1 Introduction 



In this paper we solve transport equations describing production of the quark- 
gluon plasma in strong color fields. A novel feature of our approach is the 
implementation of the boost-invariant tunneling along the hyperbolas of con- 
stant invariant time r = y/t'^ — z^ into the framework of the relativistic kinetic 
theory. In this way we generalize the previous results of Refs. [jl],|^,^,^. 

In the standard WKB description of the tunneling process P,^, the particles 
tunnel at fixed time t. They emerge from the vacuum at a certain distance 
from each other, and their longitudinal momenta vanish (in the center-of- 
mass frame of a pair). In Refs. |]l],|^,m|,0] , the WKB results were used to 
fix the longitudinal-momentum dependence of the production rates of quarks 
and gluons. The effect of the finite distance appearing between the tunneling 
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particles is more difficult to include. One assumes usually that this distance is 
small, and the non-local character of the tunneling process is neglected. This 
approximation is not always valid. For large transverse mass, the distance 
between the tunneling particles may be quite substantial and should be taken 
into account. 

The main difficulty connected with the non-local features of the tunneling 
concerns the causal properties of sequential decays. In this case one cannot 
define which pair is produced earlier or later, which leads to ambiguity in the 
determination of the decay probabilities. Nevertheless, at very high energies 
a solution to this problem exists [Q . Assuming that the pairs are produced at 
fixed invariant time r = \/t'^ — z^^ one introduces a Lorentz-invariant sequence 
of pair production, and the requirements of causality can be naturally fulfilled. 
The tunneling of particles along the hyperbolas of constant invariant time leads 
to finite longitudinal momenta of the created particles. Thus, the production 
rates in the kinetic equations should include also this extra effect. 

The boost-invariant description of the tunneling process was the main ingre- 
dient of the simulation program used to describe the space-time evolution of 
the color-flux tubes [miO]- In particular, this program was applied in the in- 
vestigations of intermittency [ll[] and soft photon emission |T^ . In the present 
paper the boost-invariant tunneling along the hyperbolas of constant invariant 
time is included in the kinetic equations. In contrast to the simulation studies 
rT],|T2[, which were restricted to the case of the elementary color flelds, we 



now deal with stronger flelds. Similarities and differences between the present 
approach and Refs. [0,|1 are discussed in detail. 



2 Pair production in strong chromoelectric fields 

2.1 Boost-invariant tunneling 



A semi-classical boost-invariant description of pair production in chromoelec- 
tric flelds was introduced in Ref. [§]. In this approach the tunneling particles 
move along the hyperbola of constant invariant time 



T= Vt2 _ ^2_ 



If the virtual particles start their motion at z = 0, see Fig. 1, the end points 
of the tunneling trajectory may be determined from the energy-momentum 
conservation laws 



Ef = Jml + p2|| = F Zf, pfw =F{tf-T). (2) 



Here Ef and pf\\ are the energy and the longitudinal momentum of the particle 
at the space-time point where it emerges from the vacuum, m±_ is the trans- 
verse mass, and F is a constant force acting on the particle (in this Section we 
shall concentrate our discussion on the case F > only). Using the rapidity 
variable y we may write 

Ef = m±cosh.yf , pf\\ = m^sinhyf . (3) 

In the analogous way we define the quasirapidity variable t], which gives 



tf = Tcoshrif , Zf = T sinh rjf. 

Equations (D, (D and §) yield @ 

m± 



sinh yf 



2Fr' 



^/ = 2 Vf- 



(4) 



(5) 



An interesting feature of Eq. (|^) is that the kinematical quantities character- 
izing the tunneling particles depend on r, i.e., pairs with different momenta 
are created at different invariant times. This is an effect of the boundary con- 
ditions which describe expansion of the system. For very large r, when the 
boundary of the field is far away from the center of the system we recover the 
standard results: pf\\ =0, zj = m±/F and tf = r. 

Although the boundary conditions change the kinematics of the tunneling 
process, the probability of tunneling per unit volume of space-time does not 
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Fig. 1. Boost-invariant tunneling of particles along the trajectory of constant invari- 
ant time r = \fW—z^ 



change [§. Thus, we may use the well established formula for the production 
rate iiJUJll 



dN 



d^x d'^p± 4:71^ 



In 1 =1= exp 



-Km I 



(6) 



Here the minus sign is appropriate for fermions (in our case for quarks and 
antiquarks) and the plus sign should be used for bosons (in our case for gluons) . 
Eq. (^ gives the rate integrated over the longitudinal momentum. The in- 
dependence may be taken into account by the following modification of (0) 



dN 



dN 



P 



d^x d^p An^ 



In 1 =1= exp 



vrrri I 



^ ( 1/ - ^ + §^/ ) • (7) 



The production rate (0) is boost invariant and reduces to the previous for- 
mula if divided by p° and integrated over py. Moreover, the rapidity of the 
particles which are produced at r/ = r/j is simply yf = ^rif, the result required 
by the condition of the boost-invariant tunneling discussed above. In the limit 
of large invariant times or large color fields we find rjf — > 0, hence the stan- 
dard formula for tunneling is recovered: the longitudinal momenta of particles 
tunneling at 2; = are zero. 



2.2 Boost-invariant variables w and v 



In the next Sections we shall use the boost-invariant variables introduced in 



Refs. P|T1 



7i = r^ = t^ — z^, w = tpw — zE, p_|_. 



and also 



V = Et — p\\ z = Jw"^ + m\u. 



(8) 



(9) 



From these two equations one can easily find the energy and the longitudinal 
momentum of a particle 



E = p' 



Q vt + wz 



P\\ 



wt + vz 



u u 

The invariant measure in the momentum space is 

dP = d-p/^ = d-pj^. 
p" V 

In addition we have 

w = Tm± sinh {y — rj) , v = Tm± cosh {y — rf) . 



(10) 



;ii) 



'12) 



Equation (|T2D allows us to rewrite the production rate in the form 



dN 
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dT "^ d^x d^p 47r3 



In 1 =F exp 




S (w — Wo) V, 



(13) 



where 



Wq = Tirtx sinh (yj ~ Vf) = — rm_|_ sinh 



Vf 



m. 



2F 



(14) 



2.3 Finite-size corrections 



Comparing Eqs. (§) and (^^ we find that the longitudinal momenta of the 
tunneling particles at 2; = are —m']_/{2FT). Thus, for large invariant times 
p\\ —>■ 0. On the other hand, in the limit of small invariant times p\\ becomes 
infinite. As we shall see in more detail in Section 4, this divergence leads to 
infinite values of color currents at r = 0, and to infinite decay rate of the initial 
chromoelectric field. The physical origin of this singularity is the possibility 
of the creation of particles at infinite values of t and z, which correspond to 
small values of the invariant time r. In practice, we always deal with finite 
systems and such tunneling cannot take place. In order to take into account 
such finite-size corrections we impose an additional condition on the tunneling 
process, namely |]10i 



2r]f = 4Arcsinh ( — — — j < Ar/. 



(15) 



Here At] is a parameter which determines the space-time region in quasirapid- 
ity, \r)\ < Ar], where the tunneling is possible. The tunneling particles should 
fit into this region, hence the distance between the emerging members of a 
pair should be smaller than Ar]. As a consequence, our final expression for the 
production rate is Q 



dN 



47^3 



In 1 =1= exp 
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-Wo)v, 
(16) 



where 9 is the step function 

e{x) = 1 for X > 0, e{x) = for x < 0. 



(17) 



^ We neglect here the boundary effects: The tunneling processes can start from a 
place close to the edge of the system and may have not enough space to fit into the 
allowed region. Such situations should be eliminated by introducing an additional 
constraint. 



We note that condition (|T^) was also used in the simulation program P,p!0| . In 
that case, the allowed region in the quasirapidity space was determined by the 
actual size of a decaying color flux tube. In the present study, the size of At] 
is suggested by the rapidity range accessible in the ultra-relativistic heavy-ion 
collisions at SPS and RHIC. In the following, we consider three typical values: 
Ar] = 4, 6 and 8. 



3 Semi-classical kinetic equations for quark-gluon plasma 



In the abelian dominance approximation, the equations for quarks, antiquarks 
and gluons have the form P, [T7|JT8[ 



p^d, ± ge, ■ F'^V^) Gt{x,p) = ^., 



{p^d^ + 



9ri.,-F"'P'^9^)G^Ax,p) = ^, (19) 



where Gf{x,p) , G^{x,p) and Gij{x,p) are the phase-space densities of quarks, 
antiquarks and gluons, respectively. Here g is the strong coupling constant, and 
i,j = (1,2,3) are color indices. The terms on the left-hand-side describe the 
free motion of the particles and the interaction of the particles with the mean 
field Ffj,iy. The terms on the right-hand-side describe production of quarks and 
gluons due to the decay of the field. The distribution functions Of and Gij 
include the spin degeneracy factors. In the numerical calculations we neglect 
the quark masses and assume that Eq. (|T8|) holds for Nf=3 flavors. We note 
that Eqs. ([I8|) and (|19D do not include the thermalization effects. The latter can 
be taken into account in the relaxation time approximation (see for example 



Refs. 19,20,2| 



The only non-zero components of the tensor F^j, = {F^^,,Fl^^) are those cor- 
responding to the chromoelectric field £, which may be written as 

£: = F30=-2— = -. (20) 

du T dr 

Here h is a function of the variable u = t"^ only (note that £ is invariant under 
Lorentz boosts along the 2;-axis). The quarks couple to the chromoelectric field 



£ through the charges p2 




(21) 

The gluons couple to £ through the charges rj^j defined by relation 

mj = e. - e,. (22) 



According to our discussion from the previous Section, the production rates 
of quarks and antiquarks in the chromoelectric field S are 



dm 



dT '"^^''^^ 
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where we have defined 
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(23) 



, (24) 



Aj = {g \ei ■ S\ - Gg) 9 {g je^ ■ S\-a, 



<ij ' 



and 



Wi 



•|^sign(e..£:). 



(25) 



(26) 



The quantity Aj describes the effective force acting on the tunnehng quarks. 
The effect of the screening of the initial field by the tunneling particles is taken 
into account by the subtraction of the " elementary force" characterized by the 
string tension aq. 



Similarly, for gluons we have 



dNi. 



dV 



'JZij{T,p^)6{w -Wij)v, 



(27) 



where 



'^'Ar,P.)-ti 


In 1 + exp — ^ 
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and 



Wi 



^ sign (rj^. ■ S) . 



2A, 



,(28) 



(29) 



(30) 



We note that the string tension of a tube spanned by gluons is three times 
stronger than that of a quark tube, a^ = Sa^ 0. 

The implementation of the boost invariance in Eqs. (|T8|) and (^) leads us to 
the following form of the transport equations H 



-^ =F g^i ■ -j^-g— = ^i^-i I^'P-l) o[wT Wi{r,p±)) , (31) 

dG ■ ■ dh dG 

Their formal solution is 

T 

Gf {r,w,p^) = Jdr' r'n, {r',p^)6{Ah, {r,r')±w- w,{t',p^)) , (33) 





Gij (r, w,p±) = / dr' r TZij {t',p±) 6 {Ahij (r, t') + w - Wij{T',p±)) , (34) 



where 

Ahiir,r') = ge,-[hiT)-hir')], (35) 

AK, (r, r') ^ grj,^ ■ [h (r) - h (r')] • (36) 



One may notice that the distribution functions (^) and ( p4D satisfy the fol- 
lowing symmetry relations 

Cr (r, w, p_l) = Gf (r, -?i;, p_l) , Gy (r, w, p_l) = Gji (r, -w, p_l) • (37) 

We note also that the time integrals in (RBf) and (53) reveal the non-Markovian 



character of the particle production mechanism: the behavior of the system 
at a time r is determined by the whole evolution of the system in the time 
interval < r' < r. 



4 Color currents 



Equations (^) and (^) show how the particles behave under the influence 
of the field. In order to obtain a self-consistent set of equations we should 
have also the dynamic equation for the field. It can be written in the following 
Maxwell form 

a^F^^(x)=r(x)+j^(x), (38) 

where ]" is the conductive current (related to the simple fact that particles 
carry color charges ej and ry^) and j^ is the displacement current (induced by 
the tunneling of quarks and gluons from the vacuum) . 



4-1 Conductive current 



The form of the conductive current is standard 

3 



r{x) = g / dPp 



NfY.^i [Gtix,p) -Gi {x,p)) + Y. nijGijix,p) 



1=1 



(39) 



Substituting the quark and gluon distribution functions (|33D and (p^ ) into 
Eq. (0), and using Eqs. ( [T0| ) and (^) we find that j^(x) has the following 
space-time structure 



j'^(^) 



'fix), 0,0, fix)] =[z,0,0,t]JiT) 



(40) 



where 



Jir) = ^ I d'p^ 
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(41) 



The use of the explicit form of the distribution functions Gi and Gij in Eq. 

n\) gives 



J{t) 
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*>J 






2 , ^2 



[A/iij- - w;ij(p_L, ^')] + P\u 



(42) 



^.^ Displacement current 



The structure of the displacement current is less obvious. One can find the 
form of j^ through the analysis of the energy and momentum conservation 
laws for both the matter (quarks and gluons) and the field 



The V = component of this equation gives for the field part 



(43) 



(44) 



and for the matter part 



^.^matter = -9^"' ' I dP p' Pu 31 



3 3 



+ fdPp^ 



N,±('-^.m.±'"'-' 



dT dT 



dT 



(45) 



where we used Eqs. (|T8D and (|1^). Using the last resuhs we may write 



p30 _ ^^ _ T730 



dt 



!■ 



F-'" -g dPp 



F^" • / dP 



3 3 



NfE 



p^Ei (dNt dN: 



^^, e, ■ F30 \ dT ■ dT 



3 p'^rjij (dNij dN,i 



E 



+ 



; 77,, ■ F30 \ dT dT 



i>3 'tj 



(46) 



Similarly, we may analyze the z/ = 3 component of Eq. (]43|) . The conclusion 
is that the field equations (^) represent a sufficient condition for the conser- 
vation of energy and momentum if the displacement current has the structure 



where 



j2,(x),0,0,j^(x) =[z,0,0,t]JD{T) 



(47) 



N 



f V^ ^i 



'^n{r) = -iY: 



r2 tiei-£ 



dw d'^p± 



'dNt dm 
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dT dT 



1 3 



+ 2 2^ „ C 



%• [ a..,a2^ (dNi, , dNj, 



dw d p± 



+ 



dT dT J 
Integration over w in Eq. (^) can be easily done 

n AT 3 „ 

JD{r) = ^Y.^ d'p^^wf+plu 7^. (r,px) 

2 -sr^ llij f j2^ 



(4J 



r2 ^ rj ■ £ 



c? PL\/wfj + p\u TZij {t,p_l) . 



(49) 



We note that the form of Eqs. (|40|) and (^Tf ) implies that both the conductive 
and displacement currents are conserved separately 



du^ix) = 0, 



duUx)=0. 



(50) 



We also note that the integral over p_j_ is restricted by condition (pISj), implic- 
itly included in TZi and 'JZij. If we did not use our finite-size correction, the 
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displacement current would diverge at small r, J' oij) ~ t~^, and the field 
equation would be singular, compare Eq. (|5l|) . In the standard case (with zero 
longitudinal momenta of the tunneling particles and with no finite-size cor- 
rections), the quantities Wi and Wij are zero, hence J" oij) ~ t~^ for small r, 
and the field equation is regular in the limit r ^ 



4-3 Field equations 



With the all substitutions required by the boost invariance, the field equation 
p8| ) may be written as 

^ = -r [J (r) + Jn (r)] (51) 

or 



dr"^ T dr 

This is an integro-differential equation for the function h(r), because the 
conductive current ^ (r) depends not only on h (r) but also on all the values 
of h(r') for < r' < r. Eq. ([5^) has to be solved numerically step by step for 
given initial values. These are taken in the form 0,^ 

h(0)=0, -^iO) = -eo = -\^,kq. (53) 

T dr V Trr^ 

Here the Gauss law has been used to determine the initial strength of the 
chromoelectric field £o in terms of the transverse radius of the color-fiux-tube 
(7rr^= 1 fm^), the string tension {ag = Sag = 3 GeV/fm), and the number of 
color charges k. Since the exchange of color charges at the initial stage of a 
heavy- ion collision leads to the color fields spanned by gluons p3[ , we assume 
that q is one of the gluon color charges rj^j. In practice we take q = 77 ^^2) 
so only the third component of the chromoelectric field, £ , is present in 
our numerical calculations. The solution of Eq. (|52D is independent of the 
initial condition for h (r) because of the cancellations connected with the 
gauge transformation which leaves £ unchanged. 

In Fig. ^ we plot the time dependence of the chromoelectric field, as calculated 
from Eqs. (^2]) and (|53D for different values of the parameter k. The maximal 
allowed quasirapidity interval Ar/ is 6 in this case. We observe the field os- 
cillations with the frequency growing with k. The shape of the oscillations is 
almost identical to those found before in Ref. [Q]. Clearly, the modification of 
the tunneling process does not affect the field behavior in this case. In Fig. ^ 
we show the time dependence of the field for the fixed initial strength, k = 3, 
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and for different values of Arj. In the three considered cases, Ar^ = 4, 6 and 
8, the field oscillations are very similar. Looking in more detail, we observe 
that with increasing A?7, the chromoelectric field decreases faster at the very 
initial stage of the process, i.e., for < r < 0.1 fm. Later the decrease of the 
field is weaker and for r ~ 0.3 fm the values of the chromoelectric fields are 
approximately the same for different values of Arj. For longer times one may 
notice that the period of the field oscillations is slightly longer for larger values 
of A?7 - a faster initial decay of the field causes a faster back reaction of the 
induced currents, and the subsequent slow down of the decay. 
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Fig. 2. Time dependence of the chromoelectric field for different values of k. The 
maximal quasirapidity interval allowed for tunneling is fixed, Arj=6. 
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5 Energy density and pressure of the plasma 



The energy-momentum tensor of the quark-gluon plasma has a structure 



rplJkV 

matter 



[e (r) + P (r)] u^'u'' - P (r) g 



flU 



where in our two-dimensional model 



«^ = - (t, z) . 

T 



(54) 



(55) 



It is important to emphasize that in our case the standard form of the energy- 
momentum tensor, Eq. (0), does not follow from the assumption of the lo- 
cal thermodynamic equilibrium, but it is a direct consequence of the boost- 
invariance. Eq. (Q) can be derived directly from the definition of the energy- 
momentum tensor with the help of the symmetry relations (^). The conser- 
vation law (p3D together with Eq. (|54D imply 



d 
dr 



'e{r) + l-£\r) 



e{T 



P(t) 



(56) 



In the absence of the fields Eq. ( |5B| ) is reduced to the Bjorken equation de- 
scribing the evolution of the energy density in a boost-invariant hydrodynamic 
model [21|. On the other hand, neglecting the expansion effects, described by 
the term on the right-hand-side of Eq. (|56|), we obtain the simple conservation 
law for the total energy of the field and matter. 

The non-equilibrium energy density e and the non-equilibrium pressure P are 



e T 



and 



P(r) 



— / d'^p± dwv 
u J 
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/ E {Gt + Gr) + 5: G., 




^ f ^2 ^ W 
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U J V 


2 


NfJ:{Gt + Gr)+J2C 





(57) 



(5J 



The time dependence of the energy density e (r), following from Eq. (0), is 
shown in Figs. § and ^. In Fig. ^ we fix the maximal quasirapidity interval 
A?7 = 6 and show the results for different strengths of the initial field, k = 2, 
3 and 5. We observe that the energy density grows very rapidly and reaches 
maximum at a fraction of a fermi. Later the energy density decreases, which 
is an effect connected with the longitudinal expansion of the system, imposed 
by the boost invariance. The maximal values of the energy density are very 
close to those found in Ref. 0. In Fig. |^ we show the results for the fixed 
initial value of the field, k = 3, and for different quasirapidity intervals, A77 = 
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Fig. 5. Time dependence of the energy density of the plasma for different values of 
Ar/ and A: = 3. 

4, 6 and 8. It is interesting to observe that with increasing values of At] the 
maximal energy density gets smaller. This behavior is connected with the time 
dependence of the chromoelectric field. For larger values of Ar/ the decay of 
the field is slower (except for the very beginning of the decay process) and the 
growth of the energy density is weaker. 

We have checked that our numerically evaluated functions £:(r) and P{t) obey 
Eq. (^61). The time-dependence of the non-equilibrium pressure P{t) is de- 
picted in Figs. 1^ and |^ One can notice that the minima of P{t) correspond 
to the extremes of the chromoelectric field. At these points the longitudinal 
momenta of the particles practically vanish and the field is the strongest. This 
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Fig. 7. Time dependence of the pressure of the plasma, k = 3 and Arj = 4, 6, 8. 

behavior reminds oscillations of a simple string. Another interesting feature 
of the process discussed here is that the ratio P (r) /e (r) oscillates around 
the mean value P/e ~ 1/3, which corresponds to the equilibrium limit. Im- 
posing condition P/e = 1/3 in Eq. ( [56| ) we may treat this expression as an 
equation determining e (the time dependent chromoelectric field enters here 
as the only input). We have found that the solutions of this equation are very 
good approximations of the exact solutions shown in Fig. | and ^. 
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6 Conclusions 



In this paper, tlie boost invariant tunneling has been incorporated into the 
framework of the kinetic theory describing production of the quark-gluon 
plasma in strong color fields. Our description of the tunneling process in- 
cludes the effect of a finite space-like distance formed between the particles. 
Such a distance appears, since the energy and momentum conservation laws 
should be satisfied locally during the tunneling. In our approach, the particles 
emerging from the vacuum have finite longitudinal momenta, which is a direct 
consequence of the tunneling along the trajectories of constant invariant time. 

We restrict the possibility of creation of pairs at very large distances (in order 
to retain the boost-invariance of the system, this condition is formulated in the 
quasirapidity space). In this way we mimic the constraints imposed by the fact 
that total energies of the realistic systems are finite. We find that the finite-size 
corrections are crucial to have non-trivial solutions of the kinetic equations. 
For the realistic finite-size corrections (as suggested by the accessible rapidity 
range in the present experiments with the ultra-relativistic heavy-ions) we find 
the solutions of the kinetic equations which are very close (qualitatively and 
quantitatively) to the results of the previous investigations: the chromoelectric 
fields oscillate, and large densities of quarks and gluons are produced in a very 
short time. The increase of the region allowed for tunneling does not lead to 
the increase of the produced maximal energy density of quarks and gluons. 
This feature of the model can be used to interpret small differences in many 
characteristics of the heavy-ion collisions at the SPS and RHIC energies. 
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